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A study of the charmless semileptonic decay B° — > 7T^£ + U£ using a large 
sample of untagged T(AS) — > BB events collected with the Belle detec- 
tor at the KEKB asymmetric e + e~ collider is presented. The branching 
fraction of the decay is obtained and from this the Cabibbo-Kobayashi- 
Maskawa matrix element \ V u b\ is determined using two different approaches. 
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1 For the Belle Collaboration 



1 Introduction 



The branching fraction for the exclusive weak decay B° — > 71 i + v provides an attrac- 
tive means for extracting the Cabibbo-Kobayashi-Maskawa (CKM) matrix element 
\V u b\. The extraction method requires knowledge of the form factor f+{q 2 ) for the 
decay, where q 2 is the square of the four-momentum transfer to the charged lepton- 
neutrino system. Traditionally this input comes from theoretical calculations such 
as those based on lattice QCD [1] [2], light cone sum rules j3] or a relativistic quark 
model [1]. The present study obtains the differential branching fraction for this decay 
in 13 regions of q 2 and from this extracts \ V u b\ using the above theoretical predictions. 
It also employs a variation on this approach which has appeared more recently [5] [6] , 
based on a simultaneous fit to data and lattice predictions, from which a value for 
\V u b\ emerges in a more model- independent way. A more extensive description of the 
study presented here is available through the arXiv e-print server [7]. 

2 Extraction of the branching fraction 

The data sample consists of 657 x 10 6 BB pairs collected by the Belle experiment [8] 
at the KEKB asymmetric e + e~ collider operating at the T(4S) resonance. B° — > 
decays are reconstructed from pairs of oppositely charged leptons (muon or 
electron) and pions, with the neutrino reconstructed from the missing 3-momentum 
in the center of mass frame (CM); p v = (|p miss |,p m i ss ), where p m i SS = —^iPi and the 
sum extends over all observed particles in the event. Signal event selection uses the 

beam-energy-constrained mass M^ c = J El eam — \p n + pe + p u \ 2 and AE = -Ebeam — 
(E n + Ei + E v ), where -EWm is the beam energy in the CM. The conditions M^ c > 
5.19 GeV/c 2 and |Ai?| < 0.1 GeV are imposed. Backgrounds are estimated using 
Monte Carlo samples of other B — > X u lv decays, B — > X c iv decays and continuum. 

A binned extended likelihood fit is used to extract the signal, in the two-dimensional 
(Mbc AE) plane and in 13 bins of q 2 . Projections of the fit are shown in Fig. [TJ After 
the q 2 distribution of the signal events is unsmeared and corrected for signal effi- 
ciency, the distribution of the partial branching fraction is calculated, with the result 
displayed in Fig. |2j Of the four models which are compared to the partial branching 
fraction distribution in Fig. [21 the ISGW2 model is excluded by the data. 

Summed over the complete range in q 2 , the branching fraction for the decay is 
determined to be B (B° -> n-£ + u e ) = (1.49 ± 0.04 (stat.) ± 0.07 (syst.)) x 10" 4 . 

3 Determination of \V u b\ 

\V u b\ can be extracted from the partial branching fraction information through the 
relation |V u6 | = \J AB(q 2 ) / (t b oAQ, where t b o = 1.525 ± 0.009 ps is the B° lifetime. 
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Figure 1: Fit projections in AE with M bc > 5.27 GeV/c 2 for low q 2 < 16 GeV 2 /c 2 
(a) and high q 2 > 16 GeV 2 /c 2 (b) regions, and corresponding projections in Mb c with 
\AE\ < 0.125 GeV in (c) and (d). Data points show statistical errors only. The 
histograms show, from top to bottom, signal B° — > tt~£ + i>, B — > X u iv^ B — > X c tv 
and continuum components. 
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HPQCD pQ 


> 16 


2.07 


± 0.57 


3.55 


± 


0.09 


± 


0.09 


+0.62 

-0.41 


FNAL [2] 


> 16 


1.83 


± 0.50 


3.78 


± 


0.10 


± 


0.10 


+0.65 
-0.43 


LCSR [3] 


< 16 


5.44 


± 1.43 


3.64 


± 


0.06 


± 


0.09 


+0.60 
-0.40 


ISGW2 [1] 


all 


9.6 


± 4.8 


3.19 


± 


0.04 


± 


0.07 


+1.32 
-0.59 



Table 1: Values for \ V u b\ extracted using different form factor predictions. The quoted 
uncertainties in \V u t,\ are statistical, systematic and arising from the uncertainty in 
AC, respectively 
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Figure 2: The partial branching fraction as a function of q 2 . Data points (closed 
circles) indicate the statistical uncertainties (inner bars) and the combined statistical 
and systematic uncertainties (outer bars). The histograms show predictions based on 
various form factor prescriptions. The curve is the result of a fit to the data using 
the BK form factor parameterization p]. 




Figure 3: The result of a simultaneous fit to the present experimental data points 
(closed circles, blue) and FNAL/MILC lattice QCD points [6 J (open circles, red) 
used to extract \V u b\. The data points have been scaled by the value of \V u b\ obtained 
from the fit. For the data points both statistical and systematic uncertainties are 
combined; for the lattice points the inner bars show statistical uncertainties and the 
outer bars combined statistical and systematic uncertainties. The functions P + and 
4> + are described in reference [6]. 
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The normalised partial decay rates A( are predicted by theory The results, which are 
dominated by theoretical uncertainties, are summarised in Tabled] for several sources 
of input for A(. 

A model- independent determination of \V u b\ is obtained by simultaneously fitting 
the branching fraction data and the MILC lattice QCD form-factor after transforming 
to the so-called "z-parameterization" [6]. The result of the fit is shown in Fig. [3], with 
the result \ V u b\ = (3.43 ± 0.33) x 1CT 3 , where the uncertainty combines both statistical 
and systematic contributions. 



4 Conclusions 

The branching fraction B (B° n-£ + u e ) = (1.49 ± 0.04 (stat.) ± 0.07 (syst.)) x 10~ 4 
has been measured and from this \V u b\ has been determined, using both model- 
dependent and model-independent approaches. The latter approach provides the best 
overall precision, producing a value \ V u b\ = (3.43 ± 0.33) x 10~ 3 where the uncertainty 
incorporates both statistical and systematic effects. 
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